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Abstract 

We perform the first complete analysis of nd scattering at next-to-next-to-leading order in chiral 
effective field theory including the corresponding three-nucleon force and extending our previous 
work, where only the two-nucleon interaction has been taken into account. The three-nucleon force 
appears first at this order in the chiral expansion and depends on two unknown parameters. These 
two parameters are determined from the triton binding energy and nd doublet scattering length. 
We find an improved description of various scattering observables in relation to the next -to-leading 
order results especially at moderate energies (i?iab = 65 MeV). It is demonstrated that the long- 
standing Ay-problem in nd elastic scattering is still not solved by the leading 3NF, although some 
visible improvement is observed. We discuss possibilities of solving this puzzle. The predicted 
binding energy for the a-particle agrees with the empirical value. 
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1 Introduction 



Effective field tlieory has become a standard tool for analyzing the chiral structure of Quantum 
Chromodynamics (QCD) at low energy, where the perturbative expansion in powers of the coupling 
constant cannot be used. The chiral symmetry of QCD is spontaneously broken and the corresponding 
Goldstone bosons can be identified with pions, if one considers the two flavor sector of the up and 
down quarks as done here. The pions are not exactly massless as it would be the case for massless u 
and d quarks, but are much ligther than all other hadrons and are therefore sometimes called Pseudo- 
Goldstone bosons. It is a general property of Goldstone bosons that their interactions become weak 
for small momenta. Chiral Perturbation Theory (CHPT) is an effective field theory which allows 
to describe the interactions of pions and between pions and matter fields (nucleons, /O-mesons, A- 
resonances, . . .) in a systematic way. This is achieved via an expansion of the scattering amplitude in 
powers of small external momenta and the pion mass. Pion loops are naturally incorporated and all 
corresponding ultraviolet divergences can be absorbed at each fixed order in the chiral expansion by 
counter terms of the most general chiral invariant Lagrangian. 

This perturbative scheme works well in the pion and pion-nucleon sector, where the interaction van- 
ishes at vanishing external momenta in the chiral limit. The situation in the purely nucleonic sector 
is somewhat different, since the interaction between nucleons is strong and remains strong even in 
the chiral limit at vanishing 3-momenta of the external nucleons. The main difficulty in the direct 
application of the standard methods of CHPT to the nucleon-nucleon (NN) system is due to the 
non-perturbative aspect of the problem. One way to deal with this difficulty has been suggested by 
Weinberg, who proposed to apply CHPT to the kernel of the corresponding integral equation for the 
scattering amplitude, which can be viewed as an effective NN potential |]. 

Following this idea Weinberg was able to demonstrate the validity of the well-established intuitive 
hierarchy of the few-nucleon forces: the two-nucleon interactions are more important than the three- 
nucleon ones, which are more important than the four-nucleon interactions and so on. 
The first quantitative realization of the above idea has been performed by Ordonez and co-workers, 
who derived the 2N potential and performed a numerical analysis of the two-nucleon system Q. To 
calculate an expression for the effective Hamiltonian for two nucleons the authors of |^] made use 
of Rayleigh-Schrodinger perturbation theory (the method is closely related to the Tamm-Dancoff 
approach @, P), which leads to a non-hermitian and energy-dependent potential. The A-degree of 
freedom has been included explicitely. The 26 free parameters, many of them being redundant due to 
the property of antisymmetry of the wave functions, have been fixed from a global fit to the low-energy 
observables. Ordonez et al. obtained qualitative fits to deuteron properties as well as quantitative fits 
to most of the scattering phase shifts up to Siab = 100 MeV. 

The property of the effective potential from Ref. of being explicitely energy-dependent makes it 
difficult to apply it to systems different from the two-nucleon one. In fact, such an energy dependence 
is not a fundamental feature of the effective interaction and can be eliminated by certain techniques, 
see e.g. Q. In ||^] we have demonstrated how to derive the energy-independent and hermitian potential 
from the chiral Lagrangian using the method of unitary transformation |^]. The advantage of this 
scheme is that it is easily extendable to processes with more than two nucleons and/or external fields. 
In we applied the above mentioned method to calculate the NN scattering observables and deuteron 
properties up to NNLO in the chiral expansion. As described in Ref. the 9 unknown low-energy 
constants (LECs) related to contact interactions and the LECs C3 and C4 related to the subleading 
vrvrNN vertices have been fixed by a fit to the Nijmegen phase shifts |11| in the ^Sq, ^Si Di, ^Pi, 
^Pq, ^Pi ^P2 F2 channels below ii^iab = 100 MeV. In contrast to Ref. we did not perform a 
global fit to the data, which due to the large dimension of the parameter space and computational 
resource limitations might not lead to the true global minimum in the x^~space and cannot easily be 
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performed. Instead we introduced an alternative set of partial-wave projected LECs and considered 
each of the above indicated channels separately having at most 3 unknown parameters in any given 
partial wave. The chiral potential at NNLO has been shown to lead to a reasonably good description 
of the NN phase shifts up to E'lab ~ 200 MeV as well as of the deuteron properties. Further we 
demonstrated that including the subleading two-pion exchange at NNLO allows to improve strongly 
the NLO results without introducing additional free parameters associated with short-range contact 
interactions, which is a good indication of consistency and convergence of the chiral expansion. For 
our choice of the LECs 01^3^4 related to the subleading ttttNN interactions see |1C]. 
The few-nucleon interactions in chiral effective field theory have been first discussed qualitatively by 
Weinberg j^. The corresponding expressions have been derived later by van Kolck, who demonstrated 
that the leading contribution, which appears at NLO in the chiral expansion, cancels against the 



iteration of the energy-dependent part of the corresponding NN effective potential [12|. Such a 
cancellation in case of the two-pion exchange 3NF has already been observed earlier jl^. Thus the 
first non- vanishing contribution to the 3NF appears at NNLO. Note that if the A-resonance is included 
explicitely, i.e. if the AA^ mass splitting is considered as a small quantity of the order of the pion 
mass, the non- vanishing contributions to the 3NF are shifted to the NLO. Note, however, that such a 
scheme is not strictly rooted in QCD because of the decoupling theorem [|l4| (but can be justified in 
the large A'^c-expansion, with Nc the number of colors). For the energy-independent potential derived 
with the method of unitary transformation one observes the vanishing of the NLO 3NF as well (as has 
also been pointed out in a different context e.g. in Ref. |jl^) and the first non- vanishing contributions 
appear at NNLO. 

In our work \1G\ we performed a complete analysis of the low-energy nd scattering at NLO in the 
chiral expansion with the NN potential introduced in H] and also calculated the triton and a-particle 
binding energies (BE's). Since no 3NF has to be included at this order and all parameters in the NN 
potential are fixed from the 2N system, the results for A > 2 systems are parameter-free predictions. 
We demonstrated a reasonably good description of the nd elastic scattering data at E'lab = 3 MeV and 
Eiab = 10 MeV as well as of some break-up observables at E'lab = 13 MeV while significant deviations 
from the data were found at E'lab = 65 MeV. The predicted value for the triton BE is in the range com- 
parable to the one based upon various modern phenomenological potentials, while for the a-particle 
BE somewhat larger deviations have been observed depending on the chosen cut-off value. Extending 
the analysis to NNLO requires, as already stated before, not only the appropriate modification of the 
NN interaction, but also the inclusion of the 3NF. In |1C] we presented an incomplete NNLO analysis 
of the 3N system based upon the NN interaction at NNLO and without inclusion of the 3NF.[| In this 
work we present the complete NNLO analysis of the low-energy nd scattering including the chiral 
3NF. We also predict the a-particle binding energy. This is the first time that the complete chiral 
3NF has been included in few-body calculations. Some pioneering steps in that direction based upon 
the hybrid approach have been done in Ref. [17|. 

Our manuscript is organized as follows. In Section II we discuss the structure of the chiral 3NF and 
demonstrate that it depends on two parameters. The partial-wave decomposition of the new terms in 
the 3NF is given in Appendix A. In Section III we discuss how these unknown parameters can be fixed 
from low-energy 3N data. Then we show our results for various elastic and break-up nd scattering 
observables as well as for triton and a-particle BE's in section IV. Conclusions and an outlook are 
given in Section V. 



From the point of view of an effective field theory it makes not much sense to include only the 2N force and to omit 
the 3NF contributing at the same order. In |^| we followed however the common trend in the field of few-nucleon physics 
and calculated various 3N observables based on our NNLO 2N potential for illustrative purposes. Such comparisons are 
helpful to identify the observables and kinematics most sensitive to the 3NF. 
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2 The Chiral Three— Nucleon— Force at Next— to— Next— To— Leading 
Order 



The chiral 3NF at NNLO is given by the two-pion exchange (TPE), one-pion exchange (OPE) with 
the pion emitted (or absorbed) by 2N contact interactions and 3N contact interactions, see Fig. |l|. 
All diagrams include apart from the leading vertices with Aj = one insertion of interactions with 
Aj = 1, where the chiral dimension is defined as 

Ai = di + ^ni-2. (2.1) 

Here di and Ui denote the number of derivatives (or pion mass insertions) and nucleon fields for a 
vertex of type i. This quantity has been first introduced by Weinberg and is especially useful in the 
few-nucleon sector. In the pion and pion-nucleon sectors one usually uses a different definition. 
The contribution from the first graph in Fig. |l] is given (in the 3N c.m.s.) by [|l2| (here and in what 
follows we use the usual notation for expressing the nuclear force: the quantity V^pe is an operator 
with respect to spin and isospin quantum numbers and a matrix element with respect to momentum 
quantum numbers): 

T/3NF_ V- { {ffj- qi){Sj ■ Qj) j^aP^a^lS „x 

" ^hk 2 ' + M2)(g, 2 + M2) '^^'^^^ ' ^'-'^ 

where iji = pi' — pi; pi {pi ') are initial (final) momenta of the nucleon i and 



iciMj , 2c3^ ^ 

f 2 + f 2 ■ 

J TV J TV 



■7 •! TV 



Here, qa = 1.276 is the axial-vector coupling constant, /^r = 92.4 MeV the weak pion decay constant 
and the ci_3^4 are the LECs from the chiral Lagrangian of the order A = 1 ||l^, which also enter the 
corresponding expressions for the subleading two-pion exchange in the 2N potential. The form ( p.2| ) 
can be shown to match with the low-momentum expansion of various existing phenomenological 3NFs 
provided they respect chiral symmetry. This issue is extensively discussed in p9t| . 
We will now derive the expressions for the OPE and contact parts of the 3NF, see also [ |12| |, and show 
that due to the Pauli principle only one independent OPE term and one independent pure contact 
term appear in the 3NF.0 

Let us start with the OPE contribution and discuss first the structure of the corresponding vrNNNN- 
vertex of dimension A = 1. After performing the non-relativistic reduction for the nucleon field 
(or, equivalently, after integrating out the lower components in the heavy-baryon formalism) one 
encounters three different structures in the effective Lagrangian (in the rest-frame system of the 
nucleons): 

= ai{N^N){N^aTN) ■ Vtt + a2(iV"^CTiV)(iVViV) • Vtt + a-iiN^arN) x {N^arN) ■ Vtt . (2.3) 

where tt and denote the pion and nonrelativistic nucleon fields, and Tj are Pauli spin and isospin 
matrices. The symbol ■ (x) denotes the simultaneous scalar (vector) product in the ordinary and 

^Similar observation for the purely short-range part of the 3NF has been made by Bedaque et al. po[ ], while I. Stewart 
pointed out that the two OPE terms in the expressions for the 3NF published in ^| are not independent from each 
other. Since these statements do not appear in the literature in a complete form we decided to demonstrate this explicitly 
here. 



4 



isospin-space. Note that the terms with derivatives acting on the nucleon fields are ehminated by 
partial integration. The corresponding 3N force at NNLO is of the form 



^OPE Of (*c • Sk) ^ ,0 ■ \ "1 Si'^i + "2 OiTj + as {Oi X <Tj)(Ti X T j) \ . (2.4) 



Since we treat nucleons as identical particles, the few-nucleon states are antisymmetric. For these 
antisymmetric states the operators Vqp^ can be further simplified. Because the force is symmetric 
with respect to an interchange of particles i and j, the relation 

\^) = A, \^) = v^r^ A, 1^) (2.5) 

holds and therefore one can work equally well with an antisymmetrized force. Here Aij is the anti- 
symmetrization operator in the space of two nucleons i and j, which reads: 

A, = , (2.6) 

where Pij is the corresponding permutation operator, Pij\ij) = given by 

„ 1 + Bi ■ (Tj 1 + Ti ■ Ti , , 

Pii = — . 2.7) 

In addition, one has to interchange the corresponding nucleon momenta. It is an easy exercise to 
apply the antisymmetrization operator Aij to that pair ij of the 3NF in eq. ( |2.4D which interacts via 
the contact terms and to see that all three different structures lead to the same expression. 
In the case of the purely contact 3NF without derivatives we proceed in an analogous way. The most 
general structure of such 3NF which satisfies the usual symmetry requirements (rotational and isospin 
invariance, parity invariance and invariance under time reversal transformation) is given by 



3NF _ 

ii^i+k 



^ont 



X! ] /?1 + /32 <7i • Oj + PsTi- Tj + /34 {ai ■ aj){Ti ■ Tj) + /35 {ffi ■ Bj){Tj ■ Tk) (2.^ 

+ (pi X aj] ■ ak){[Ti x Tj] • Tfc)| . 

The antisymmetrization operator Aijk in the space of three nucleons can be expressed as: 

_ (l + P,,P,fc + P,fcP,fc) (1 - Pjk) ,^ , 



Acting with the operator Aijk on the 3NF in eq. ( |2.8| ) and performing a straightforward, but somewhat 
tedious simplification one ends up with a single structure just as in the previously considered case. 
We thus have shown that it is sufficient to consider only one OPE and one pure contact term in the 
chiral 3NF at NNLO, since all other terms have due to the Pauli principle precisely the same effect 
on the S-matrix. In what follows, we will use the following form for these 3NF contributions: 

« = - E S ^ (r. ■ r,) (a. ■ q,) , (2.10) 

i^j^k^J-^ +^^^7^ 

V!^ = \i:E{rrrk). 
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where D and E are the corresponding LECs from the Lagrangian of dimension A = 1: 



C^^'^ = --^ (Aft AT) (N^arN) ■V7V--E {N^^N) (N^tN) ■ (iVVTV) , (2.11) 
Note that dimensional scaUng arguments ahow one to express the LECs D and as |2l| 

where co and ce should be numbers of order one and A-^ is the chiral symmetry breaking scale of the 
order of the p meson mass. Here and in what follows we use A-^ = 700 MeV. It has been demonstrated 
in [^] that all corresponding numbers for 2N contact interactions at NLO and NNLO are natural 
for the cut-off values considered. It should also be understood that a more precise analysis of the 
naturalness would require also taking into account symmetry factors in the Lagrangian as well as 
additional factors resulting from insertions of spin and isospin matrices.^ 

3 Fixing the Parameters of the Three— Nucleon— Force 

We now proceed to fix the unknown LECs cd and ce from 3N low-energy observables. To that aim we 
solve the 3N Faddeev equations for the bound state and for nd scattering. They have the well known 



form [g3|, |2J 



V; = Go t P V + (1 + Go t) Go V^Sf i^ + P)^, (3.13) 

in case of the bound state. Here V^^F that part of the three-nucleon force which singles out one 
particle (here particle 1) and which is symmetrical under the exchange of the other two particles. The 
complete 3NF is decomposed as 

^3NF = V^3NF + ■^"3NF + ^NF ■ (3-14) 

Further, tp denotes the corresponding Faddeev component, t is the two-body t-operator. Go = \/{E ~ 
Hq) is the free propagator of three nucleons and P is a sum of a cyclical and anticyclical permutation 
of the three particles. In case of nd scattering we follow our by now standard path |2^, ^ and firstly 
calculate a quantity T related to the 3N break-up process via the Faddeev-like equation: 

T = tPc^+{l + tGo) y^NF (1 + P) + t PGo T + (1 + t Go) V^^f (1 + P) Go T , (3.15) 

where the initial state 4> is composed of a deuteron and a momentum eigenstate of the projectile 
nucleon. The elastic nd scattering operator is then obtained as 

U = PGq^ +PT + Vi^l (1 + P) (1 + Go T) , (3.16) 
and the break-up operator via 

Uo = {l + P)T. (3.17) 
These equations are accurately solved in momentum space using a partial wave decomposition. For 



details see 27, 28]. The corresponding partial wave decomposition of the chiral 3NF is given in 



the appendix. The equations ( 3.13| ) and ( |3.15| ) have to be regularized, since the expressions for the 



Such factors can be calculated from expressions of the 3NF. For example, the antisymmetrized expression of the 



third term in eq. (2.8) is 3 times smaller than the one of the fourth term, which has two additional insertions of the Pauli 
spin matrices. 
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3NF (U) and (|2lo| ) are only meaningful for momenta below a certain scale. We regularize the y3NF 



in the way analogous to the one adopted in the analysis of the two-nucleon system [^] : 

V'^'^ip, q; p', q') - fnip, q) V'''^ {p, q; p', q') fdp', q') , (3.18) 

where p and q {p' and q') are Jacobi momenta of the two-body subsystem and spectator nucleon 
before (after) the interaction. The regulator function fR{p,q) is chosen in the form 



fR{p,q) = exp 



J 



(3.19) 



so that it coincides with the exponential function fji^°^{p) of Ref. for q = 0. Clearly, this is not 
the only possible choice for that function. The final results for low-energy observables are insensitive 
to the choice of the regulator function provided that it does not violate the appropriate symmetries. 
Note that the values of the LECs cd and ce „run" with the cut-off A to compensate the changes in the 
observables, which are cut-off independent (up to the accuracy at the order in the chiral expansion). 
The dependence of the LECs on the cut-off A is governed by renormalization group equations, as it 
is always the case in quantum field theory. We choose A in the 3NF equal to A in the NN interaction. 
The following study has been carried through with the minimal and maximal momentum cut-offs, 



A = 500 and 600 MeV, for which our NN force has been defined in |1C]. 

The low-energy constants cd and ce enter the expressions for the chiral 3NF at NNLO. The constant 
Ce can only be obtained from 3N data, while cd can be best determined in the 3N system or, for 



larger momentum transfer, in pion production in NN collisions |28]. One important part of this work 



is to outline a feasible way to fix these parameters. We will now show that the LECs cd and ce 
can be determined using the '^H BE and the nd doublet scattering length ^a„d, which are bona fide 
low-energy observables. Since for the time being we have no nn and pp forces at our disposal (these 
have been calculated in chiral EFT to NLO so far ||3^) and both observables we are interested in are 
known to depend on the difference between np and nn forces, we decided to use np/nn corrected data 
as input to our fitting procedure. To this aim, we compare results using phenomenological forces with 
the proper np and nn forces and with a np force only. Several combinations of NN and 3N forces have 
been adjusted to describe the triton BE (see [^]). We used AV18 |32| augmented by the Urbana-IX 



3NF pi and CD-Bonn 2000 pi augmented by the TM99' 3NF p|]. These models come along with 



nn forces, which are adjusted to the nn scattering length. Replacing these nn forces by the np ones, 
we find an increased binding energy of 8.65 and 8.72 MeV, respectively. From those we estimate a np 
corrected „experimental" pseudo BE of 8.68 MeVj^ 

The „experimental" pseudo value for "^Qnd has been determined using the NN force CD-Bonn alone. 
The corresponding shift of ^a„d is -0.19 fm. This together with the experimental value ^a„d = 0.64 it 
0.04 fm leads to the „experimental" pseudo value ^o„d = 0.45 ± 0.04. It should be understood that 
the uncertainty in the estimated pseudo value of the scattering length is even larger due to the error 
in the shift resulting from replacement of the nn force by the np one. We however refrain from further 
discussion of that issue. 

For the chiral interactions at NNLO two unknown LECs enter into the 3N bound state Faddeev 
equation: cd and ce- Both affect the BE strongly. Imposing the condition that the Hamiltonian 
describes the pseudo BE, we find a correlation between both LECs, which is displayed in Fig. ^. The 
unsymmetric interval shown in this figure is a consequence of the fact that the doublet scattering 
length favors positive values for cd- The correlations have a very different behavior for both cut-off 



^"in fact, it would also be sufficient to make an estimation of the isospin brealting effects based upon purely 2N forces 
at tlie level of precision of NNLO. 
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values. For A = 500 MeV the functional form turns out to be nearly linear. This is not the case for 
A = 600 MeV. Later on we will demonstrate that this different behavior of the correlations for both 
cut-off values does not show up in observables. 

One needs a second condition to fix both LECs uniquely. The nd doublet scattering length ^a^^ 



is known to be correlated with the '^H BE. This correlation is known as the Phillips line |36|. We 
investigated it in the context of chiral nuclear forces. It turns out that the scattering length depends 
on Co even if ce is chosen according to the correlation in Fig. |^ with the fixed value for the triton 
BE. This indicates that the correlation between the doublet scattering length '^and and the BE 
is not exact. In fact, already for conventional NN and 3N forces, there was a slight scatter around 
an average line correlating the and '^Und values for different nuclear force^ [|3^. The Phillips 
line has recently been rediscovered within pion-less EFT |38|,|39|. At LO and NLO in the pion-less 



EFT the 3NF is given by a single contact term without derivatives and thus depends on just one free 
parameter. The Phillips line results from variation of this parameter and is in agreement with results 
based upon phenomenological interactions. Going to higher orders in the low-momentum expansion 
one encounters contributions to the 3NF with more derivatives and the exact correlation between 
'^ttnd and BE observed at LO and NLO is broken, see |39] for more details. As discussed above. 



in the EFT with explicit pions the first nonvanishing 3NF at NNLO already depends on two free 
parameters and thus the Phillips line is already broken at this order in the chiral expansion. This 
allows to determine cd (and at the same time ce) by a fit the the "experimental" pseudo datum for 
the doublet scattering length. In Fig. |3| the grey horizontal band indicates the scattering length range 
in agreement with the experimental error bar. Our theoretical predictions for A = 500 and 600 MeV 
are shown against cd. We read off from Figs. ^, ^ the following values: 

CD = 3.6, CE= 0.37, A = 500 MeV, 

CD = 1.8, Cij = -0.11, A = 600 MeV. (3.20) 



Notice that the sign of the determined LEG cd agrees with the one found in |29] from P-wave pion 
production in the proton-proton collisions. Note that for comparing our results with the ones of [29| 
one should take into account different conventions with respect to gA- We are aware of the fact that 
we can only obtain a first estimate of cd and ce- The most important uncertainties are the errors due 
to the np force corrections and the experimental error bar of ^a^^. In principle the errors due to these 
uncertainties with respect to observables could be estimated performing calculations with several cd 



and Ce combinations consistent with these error bars. In view of upcoming new data for ^a„rf |4C, 41 



and work on the isospin breaking in our formalism, we postpone such an analysis. In summary we 
emphasize that the breakdown of the Phillips line correlation enables us to determine the LECs from 
the 3N BE and the nd doublet scattering length. The result is a parameter free 3N Hamiltonian. In 
the next sections we will investigate the results for the 4N bound state and 3N scattering based on 
this Hamiltonian. 



4 Predictions for Three— and Four— Nucleon Systems 

We start with the prediction for the a-particle BE. This is based on the solution of Yakubovsky 



equations [42| as described in 31]. The results are fully converged and accurate to 2 keV for the 
3N and 20 keV for the 4N system. The convergence with respect to partial waves is much faster for 
the chiral interactions than for the convential ones. This is a consequence of the momentum cut-offs. 



^^One should be careful by looking at results which appear in the literature and sometimes indicate quite a strong 
deviation from the Phillips line. Especially earlier calculations have often been performed with not phase-equivalent 
potentials and with restricted accuracy. 
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T 
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NLO 


500 


-8.54 


30.76 


-39.30 




-29.57 


61.42 


-91.00 






600 


-7.53 


39.24 


-46.77 




-23.87 


77.61 


-101.47 




NNLO 


500 


-8.68 


31.07 


-39.43 


-0.318 


-29.51 


61.83 


-89.59 


-1.753 




600 


-8.68 


34.44 


-42.41 


-0.712 


-29.98 


71.49 


-97.44 


-4.025 


"Expt" 




-8.68 








-29.8 ±0.1 









Table 1: '^H and ^He BE at NLO and NNLO of the chiral expansion (for the cut-off range considered 
throughout) compared to "experimental" pseudo BE (see text). Apart from the BE's E (in MeV), we 
also give the kinetic energies T (in MeV) as well as expectation values of 2N and 3N forces Vat at and 
VsAT, respectively, (in MeV). 



which suppress the high momentum components exponentially. The calculations of the binding energy 
for the chiral interactions are truncated at a two-body total angular momentum in the subsystem of 
irnax = 6 for the 3N system. For the 4N system we truncate the partial wave decomposition by 
the restriction that the sum of all three angular momentum quantum numbers is below Z™'^ = 10. 
Calculations for conventional forces require = 14 (for details see [pl|]). 

Before we comment on our results for the BE's, we need to define a Coulomb and np corrected a— 
particle BE. Again, based on AV18+Urbana-IX and CD-Bonn+TM99', we calculated BEs for the 
a-particle of 28.5 MeV and 28.4 MeV. Replacing the pp and nn forces by np forces and omitting 
the Coulomb force, the BE's change to 29.9 MeV and 30.0 MeV. From these results we estimate an 
average change of the BE of 1.5 it 0.1 MeV. The experimental a-particle BE is 28.3 MeV. Thus we 
compare our results for the chiral interaction to an "experimental" pseudo BE of 29.8 it 0.1 MeV. 
In Table |l] this value is shown together with the "experimental" pseudo BE for in comparison to 
the NLO and NNLO results for '^H and ^He. The BE is in general very sensitive to small changes 
of the interaction, as it comes out as the difference of the large kinetic and potential energies. As a 
consequence, we found a rather large dependence of the BE's on the cut-off at NLO [|l6| (~19 % for 
the a-particle). At NNLO the BE agrees with the "experimental" value by construction. Because 
of the strong correlation of 3N and 4N BEs, known as Tjon-line |44|, one can expect a rather small 
cut-off dependence of the a-particle BE, too. This is indeed the case. However, we would like to 
mention that 3NFs break this correlation |31] and we observe a cd dependence of the a-particle BE 
(1.5 MeV change in the range cd = —1.5.. + 1.5 for A = 500 MeV). We are also very encouraged by 
the fact that the a-particle BE for both cut-offs comes out close to the "experimental" value. Note 
also that no 4NF contributes at NNLO. Therefore all predictions for ^ > 3 at NNLO are parameter 
free. 

Additionally, we list the expectation values of the different parts of the Hamiltonian in Table |l[ It is 
important to realize that those quantities are not observable. We see that the relative contributions of 
the NN and 3NF parts are comparable in the 3N and 4N system. We also observe that the ratio of NN 
and 3NFs strongly depends on the cut-off chosen. This is elaborated in more detail in Table ^, where 
the 3NF expectation value is split in the contribution from the 2tt exchange (c-terms), Ivr exchange 
(D-term) and contact term (E-term). The contributions of D- and E-term cancel to a large extent in 
both nuclei and for both cut-offs. The change in sign of the E-term changing from A = 500 MeV to 
A = 600 MeV has to be expected, since the ce changes its sign, too. More surprising is the change in 
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A 




^He 


c-terms 


D-term 


E-term 


all 


c-terms 


D-term 


E-term 


all 


NNLO 


500 
600 


-0.39 
-0.73 


0.81 
-0.12 


-0.74 
0.13 


-0.32 
-0.71 


-2.00 
-3.81 


3.93 
-0.84 


-3.69 
0.63 


-1.75 
-4.03 



Table 2: Contribution of the different terms of the 3NF to the complete 3NF expectation value for 
and ^He. All energies are given in MeV. 



sign for the D-term. This has to be caused by a qualitatively different action of the D-term operator 
on the wave functions for A = 500 MeV and A = 600 MeV. 

We now switch to scattering observables. Most of the 3N scattering data have been obtained for the 
pd system. In the case of scattering the isospin breaking effects in the nuclear force are believed to 
be of minor importance. We have checked this assumption explicitly for elastic scattering observables 
using the CD Bonn potential with np and nn and with only np forces to evaluate corresponding 
effects. Only in two cases at 3 MeV, namely for T20 (at forward angles data are shifted upwards) 
and for T21 (data are shifted downwards at angles below 120°), significant effects were found. The 
np-force corrections are small for all considered elastic scattering observables at 10 MeV and nearly 
invisible at 65 MeV. Therefore we refrain from correcting data for this effect. In contrast, there are 
visible Coulomb corrections necessary at these energies. We are not able to take the Coulomb force 
into account in the 3N continuum. For the Coulomb corrections we rely on the work of the Pisa 
collaboration, who can calculate low energy scattering observables based on the full AV18 interaction 
including the Coulomb force |^, ^] . The difference of these full calculations and calculations without 
Coulomb force serves as our estimate of the Coulomb corrections. In the following, all pd elastic 
scattering data at 3 MeV and 10 MeV have been corrected by this amount. For 65 MeV we did 
not correct the data assuming that Coulomb corrections are small except in forward direction. For 
the break-up we refrained from any corrections because of the lack of reliable theoretical calculations 
taking the Coulomb force into account. 

The nd scattering observables have been studied very intensively using the modern phenomenological 
interactions ||2^, 46]. In general, the description of the data by these models is very good at low 
energies with a few prominent exceptions. The most famous one entered the literature as Ay puzzle 



[47, 48 1, which is related to the fact that this observable is underpredicted in the maximum by realistic 



high-precision models of the nucleon interactions. In this paper we do not compare the new results to 
traditional ones. That has been done in |1C] for the NLO and NNLO interactions without 3NF part. 



In this case, the NNLO interaction compares quite well to the results |49| based on the highly accurate 
phenomenological forces. In the following we would like to concentrate on the complete analysis at 
NNLO. 

In Figs. ^, ^ and ^ we show a comparison for few selected elastic scattering observables at 3, 10 and 
65 MeV, respectively. The left column shows our results for NLO |5C] in comparison to the data and 
the right column the new NNLO results compared to the same data. The bands are given by the 
cut-off variation in the range from A = 500 and A = 600 MeV. They may serve as an estimation of 
the effects due to neglected higher order terms in the chiral expansion. 

The differential cross section is presented in the first row of Figs. ^ and |6|. Additionally, we give 
a more detailed look at the cross section minima in Fig. ^. At 3 MeV and 10 MeV we see that 
NLO and NNLO predictions overlap. The cut-off dependence is already small at NLO and nearly 
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vanishes at NNLO. This strong reduction of the cut-off dependence of this observable at NNLO is 
expected and can easily be understood. Indeed, at least at low energy the differential cross section is 
dominated by the 2N S-waves. The situation is more interesting at 65 MeV. In the minimum of the 
cross section one observes large differences between the NLO and NNLO results (also to the incomplete 
NNLO calculation, see The cut-off dependence of the NLO results is more visible than at lower 

energies, and is again strongly reduced at NNLO. The NNLO results are in agreement with the data 
except for forward directions, which are sensitive to the Coulomb force. Note that the improvement 
at NNLO is not only due to the fact that the NNLO 2N potential leads to a much more accurate 
description of the data especially at moderate energies [^], but also due to the chiral 3NF. This is 
demonstrated by the dotted line in the lower panel of figure 0, which corresponds to cd = —3.0 at 
A = 500 MeV (and a ce chosen appropriately to reproduce the 3N binding energy) .p^ For this value of 
CD the prediction in the minimum is in disagreement with the data. It is gratifying to see that fixing 
the LEG CD from the scattering data at zero energy we are able to describe the cross section minimum 
at 65 MeV. We consider this to be an important indication of consistency in the determination of the 
LECs CD and ce- 

As already pointed out before, the most problematic observable of nd elastic scattering is Ay, which 
is shown in the second row of Figs. |, I and |. First of all we would like to stress that vector and 
tensor analysing powers are defined as differences of polarized cross sections and are rather small at 
low energies, so that larger theoretical errors for these observables have to be expected. At energies 
3 and 10 MeV we see visible deviations of our predictions for Ay from the data for both NLO and 
NNLO. It is well known |51] that this observable is extremely sensitive to the ^P-wave phase shifts 
in the NN system. Although at NLO chiral predictions at 3 MeV seem to be in agreement with the 
data, this cannot be considered as a final result in chiral EFT. Indeed, the ^P-wave phase shifts in the 
NN system are only described at low energies with an accuracy of about 5% |1C], which indicates that 
large corrections to nd Ay at higher orders in the chiral expansion are possible. At NNLO the ^P-wave 
phase shifts come out with a significantly smaller error of about 2% (at Tiab = 10 MeV) and we found 
in [10 1 that Ay is underpredicted if one only uses the 2N forces, just as in the case of high-precision 
potential models. As one sees from Figs. |^ and |5|, we do not solve the Ay puzzle performing the 
complete NNLO analysis and including the 3NF. It is important to stress that in principle, one could 
try to solve this puzzle by the NNLO chiral 3NF. Indeed, instead of fixing the unknown LECs cd and 
C£; to the triton binding energy and the doublet scattering length, one could think about requiring a 
good description of Ay at, say, 3 MeV as being one of the two conditions needed. We found, however, 
that Ay is not very sensitive to the choice of the LECs cd and ce, if the two are adjusted to reproduce 
the triton binding energy. We were not able to find values of these coefficients in the natural range, 
which would simultaneously describe Ay at 3 MeV and the triton BEJ^ This indicates that higher 
order effects are probably still important for this observable. On the other hand it can be also a hint 
that the Ay puzzle is related to an insufficient knowledge of the low energy ^Pj NN phases [52|. It 
will be interesting, as a next step, to include l/rriN corrections to the interactions and to study their 
effect especially on Ay. At 65 MeV the decription of Ay is much better, which is in agreement with 
results based upon the phenomenological nuclear forces. 

The lower three rows of Figs. |^, ^ and ^ show the tensor analyzing powers T20, T21 and T22. At 
3 MeV and 10 MeV, we find in general that the NNLO predictions stay within the NLO band. The 
cut-off dependence clearly shrinks, which is a good indication of convergence of the chiral expansion. 

^^This value of cd is exclude d by the observed value of the doublet scattering length. 

^^Huber et al. pointed out in pj], that Ay is sensitive to the choice of cd- This is not in contradiction with the above 
discussion. We also observed a sensitivity to independent variations in co and Ce, which is, however, strongly reduced 
as soon as one requires to reproduce the triton BE. 
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Alltogether the agreement with the data is good except for T20 and T21 at 10 MeV.|^ Notice that similar 
results have been reported in |53] based upon the combination of the AV18 2N and the Urbana IX 
3N forces, where these observables have been calculated in the pd system and the Coulomb force has 
been taken into account. At 65 MeV, the situation is comparable to the one for the cross section. 
While the NLO predictions at this energy deviate significantly from the data, the NNLO results are in 
a much better agreement. Unfortunately, the quality of the data does not allow to draw more precise 
conclusions and especially here new high-precision data are needed. 

Let us now switch to break-up observables. We performed calculations at two energies, 13 and 
65 MeV, where a lot of pd data exist. As already pointed out above, there are no reliable Coulomb 
corrections available for the break up. Therefore we show the non-corrected pd data in comparison 
to our nd calculations. Note that for the space star configuration at 13 MeV presented in Fig. ^, it 
is shown that the nd and pd cross section data strongly deviate indicating that Coulomb effects can 
become important at least in some configurations. Presumably, the Coulomb corrections are smaller 
at 65 MeV. 

At 13 MeV we demonstrate in Fig. ^ chiral predictions for the cross section in the often investigated 
final-state interaction peak, quasi-free scattering and space-star configurations, which have also been 
considered in the NLO analysis [|l6|. For a general discussion on various break up observables and 
configurations the reader is referredd to Ref. |25|. As demonstrated in Fig. |^, the NLO and NNLO 
results essentially agree at 13 MeV. They describe the configuration dominated by FSI peaks quite 
well (for a more elaborated procedure the angular openings of the detectors have to be taken into 
account, see [25|). The present theory for the break-up configuration including a QFS geometry fails 
in the central maximum. The reason might be Coulomb force effects. The third configuration, the 
so called space-star, is one of the long standing puzzles of 3N scattering |54, 55, 56, 57, 58]. Similar 
to the phenomenological interactions, we even fail to describe the nd data, pd and nd data are quite 
different and it remains open whether Coulomb corrected data would fall on the theory. We observe 
[25 1 the tendency that with conventional NN forces theory is already rather close to the pd data at 
19 MeV and even closer at 65 MeV. This suggests that presumably the discrepancy to pd data is due 
to Coulomb force effects. 

At 65 MeV we decided to present the results for the same configurations as the ones studied recently in 

in order to enable a direct comparison between 
and include, in addition to the cross 



the context of phenomenological nuclear forces [54, 
these two different approaches. We follow the lines of |54, 
section, also Ay. The situation at 65 MeV seems in general to be very promising as documented in 
Figs. ^20. The improvements in the description of the data in going from NLO to NNLO are quite 
impressive. It is interesting that sometimes in case of Ay the band width at NNLO is still relatively 
large, which indicates that this observable might get significant corrections at higher orders. In Figs. |TT| 
and 16 we fail to describe the cross section in part of the S'-range. The reason is not known to us. 
These observables also change visibly, when going from NLO to NNLO. Here we cannot claim that 
convergence with respect to the chiral expansion is reached at NNLO. For one of the configurations 
(see Fig. |l^, the step for Ay going from NLO to NNLO is dramatic and better data would be very 
welcome. Finally we point to two more cases in Figs. 13 and 14, where the band width in the cross 
section shrinks nicely going to NNLO and where the agreement with the data is quite good. 
In view of the quite good description of the Nd elastic and break up data at 65 MeV at NNLO and 
of the good description of the NN data up to 200 MeV, we are optimistic and expect to be able to 
describe the data at NNLO in the energy regime towards 100 MeV. From investigations based on 
phenomenological interactions |4£, 54], we expect that there 3NF effects become clearly visible at 
these higher energies. In addition, observables at these energies will probably be more sensitive to the 

^^Significant deviation from the data in the minimum of T21 can be observed at 3 MeV as well, if np-force corrections 
are taken into account. 
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structure of the 3N interaction. For example, the sensitivity of the cross section minimum to the value 
of CD observed at 65 MeV is expected to be magnified at higher energies. Therefore we call for more 
data at intermediate energies, which could be compared to predictions of chiral EFT. Notice also that 
the existence of such high-quality data will be of a crucial importance for higher order calculations, 
where more parameters in the 3NF will have to be fixed from the data and a better accuracy in the 
theory will be reached. 



5 Summary and Conclusion 

In summary we applied for the first time the complete chiral EFT interaction at NNLO to the 3N and 
4N bound states and to 3N scattering. We reexamined the 3NF of the chiral interaction at NNLO and 
used antisymmetrization to eliminate all parameters except two. We showed that these two parameters 
can be determined from the BE and the ^a„rf scattering length. For the time being the accuracy 
of the scattering length is not sufficient to perform a precise determination of these two parameters. 
However, the favorable description of nd scattering data indicates that the values chosen in this work 
are in a reasonable range. 

We showed that the obtained parameter free Hamiltonian leads to a good description of the a-particle 
BE. The theory thus seems to be applicable to this densely bound system. It will be interesting to 
apply the 3N Hamiltonian to light nuclei, e.g. using the no-core shell model technique [^, |60[ . 
Overall we observe a good description of the data at NNLO. Most of the low energy elastic scattering 
data (at 3 and 10 MeV) are described at both orders NLO and NNLO showing convergence of the 
chiral expansion and agreement with the data. Ay turns out to be a problematic observable as there 
is still no agreement with the data and the predictions for NLO and NNLO disagree. Whether this 
will be cured by 1/m corrections has to be studied in a forthcoming paper. 

At 65 MeV the situation is also very promising. In general, we observe that the NNLO predictions 



move towards or onto the data, while the NLO results deviate significantly from the data. In Ref. [10| 
we found that the NNLO interaction can describe the NN phase shifts up to energies of 200 MeV 
neutron lab energy. Here we see that the extension of the energy range going to NNLO for the 
two-body system is continued in the few-body systems. 

This study was based on the systematic expansion of the nuclear force according to chiral perturbation 
theory applied to the NN potential. We emphasize that the favorable agreement with the data, the 
stability of our predictions when going from NLO to NNLO, observed in most cases, as well as the 
decreased cut-off dependence of the NNLO results indicate consistency of our calculations. New nd 
data in the energy range between 65 MeV and 100 MeV are highly welcome and would allow to draw 
quantitative conclusions on the range of validity of the NNLO approximation as well as to probe the 
spin structure of the leading 3NF. Such data would also be of a crucial importance for extending the 
analysis to higher orders. 

In the next steps, we have to take into account the isospin breaking of the nuclear force. Together 
with upcoming new data for the doublet nd scattering length a much more accurate determination of 
the 3NF parameters will than be possible. 
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Appendix A. Partial wave decomposition of the chiral 3NF 

Here, we give the explicit formula for the partial wave decomposition of the chiral 3NF. Since the 



partial wave decomposition of the TPE 3NF is already discussed e.g. in |28|, we only concentrate here 
on the remaining contributions to the NNLO 3NF due to the OPE and contact term in eq. ( |2.10| ). 
For general details on the partial wave decomposition in the 3N system the reader is addressed to 



Ref. 1 22]. As already pointed out before, we usually decompose 3NFs into three parts according to 
eq. (13.141) . In the following we give expressions for one such part K^p^p. For the OPE term we find: 

9DgA, 



— -2-(47r)^(5jj/ 6mM' ^TT'^MTM^^lO^l'O^sjSs'j' 1 + (-l)"*^"* 





J2 ykiillki, 000) 

fci=0,2 



(A-1) 



1 



V(2/i)!(2/2)! 
(A;/iA', 0)(/cZ2A, 0), 



where 



1 q2 



(A-2) 



Here Q ^ y + q' — 2qq'x and Pk{x) is a Legendre polynomial. Further, p and q {p' and q") are 
relative initial (final) Jacobi momenta in the pair jk, j, k ^ i, and of the nucleon i with respect to the 
pair jk, respectively. / and A (/' and A') denote the initial (final) relative orbital angular momenta 
within the pair jk, j, k ^ i, and between the nucleon i with respect to the pair jk. The initial (final) 
spin of the subsystem jk, j, k ^ i,\s denoted by s {s'). In addition, / and s {I' and s') are coupled to 
the total subsystem angular momentum j [j'), and A (A') and Sj = ^ to the total spectator angular 
momentum I (/'), which finally combine to J {J') accompanied by M {M'). The total isospin quantum 
numbers T Mt {T' M^) are constructed analogously: \{t^)TMT) '^'-^'^ 
a convenient abbreviation 

f = 2/ + 1 . 



(l(t'i)r'M^)). We also introduced 



(A-3) 



For the contact term in the second line of eq. ( 2.10| ) we find: 



i(Ma|^3NF, cont\p' <l'^')i 
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FIGURES 




1 2 3 



Figure 1: Three-nucleon force at NNLO: TPE, OPE and contact interaction. Solid 
and dashed lines are nucleons and pions, respectively. Heavy dots denote leading 
vertices with A, = and solid rectangles correspond to vertices of dimension Aj = 1. 
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Figure 2: Correlation between the LECs ce and cd after adjustment to the triton 
pseudo BE. 
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Figure 3: nd doublet scattering length and as function of the constant cd 
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Figure 4: nd elastic scattering observables at 3 MeV at NLO (left column) and 
NNLO (right column). The filled circles are nd pseudo data based on ^ while 
the filled triangles are true nd data |63]. The bands correspond to the cut-off 
variation between 500 and 600 MeV. The unit of the cross section is mb/sr. 
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Figure 5: nd elastic scattering observables at 10 MeV at NLO (left column) and 
NNLO (right column). The filled circles are nd pseudo data based on |p3, |6^, |63] 



while the filled triangles are true nd data |66]. The bands correspond to the cut-off 
variation between 500 and 600 MeV. The unit of the cross section is mb/sr. 
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Figure 6: nd elastic scattering observables at 65 MeV at NLO (left column) and 



NNLO (right column). The filled circles are pd data |61, 67]. The bands correspond 
to the cut-off variation between 500 and 600 MeV. The unit of the cross section is 
mb/sr. 
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Figure 7: Minima of the cross section (in mb/sr) of elastic nd scattering at 3 MeV 
(upper panel), 10 MeV (panel in the middle) and 65 MeV (lower panel) at NLO and 
NNLO. The filled circles are nd pseudo data at 3 and 10 MeV and true pd data at 
65 MeV. The bands correspond to the cut-off variation between 500 and 600 MeV. 
The dotted line at 65 MeV shows the NNLO result with cd = -3.0 and A = 500 
MeV. 
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Figure 8: nd break up cross section in [mb MeV^"*^ sr^^] along the kinematical locus 
S (in MeV) at 13 MeV in comparison to predictions at NLO (light shaded band) and 
NNLO (dark shaded band) in chiral effective field theory. In the upper row a final 
state interaction configuration is shown, in the middle one a quasi-free scattering 
configuration (both in comparison to pd data) and in the lower one a space star 
configuration (upper data nd, lower data pd). The precise kinematical description 
can be found in Ref. |2^]. pd data are from []68| , nd data from |]69| , [fo| . 
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Figure 9: nd break up cross section in [mb MeV~^ sr~^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 



field theory. Symmetric space star configuration is shown, pd data are from [71 
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Figure 10: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 



field theory. Symmetric forward star configuration is shown. pd data are from |71| 
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Figure 11: nd break up cross section in [mb MeV"^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 



field theory. Quasi-free scattering configuration is shown, pd data are from |71| 
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Figure 12: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 



field theory. Quasi-free scattering configuration is shown, pd data are from |71| 
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Figure 13: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Collinear configuration is shown, pd data are from ||7^ 
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Figure 14: nd break up cross section in [mb MeV"^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Collinear configuration is shown, pd data are from ||7^ 



30 





S [MeV] 



Figure 15: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Collinear configuration is shown, pd data are from l?^ 
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Figure 16: nd break up cross section in [mb MeV"^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Collinear configuration is shown, pd data are from l?^ 
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Figure 17: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Unspecific configuration is shown, pd data are from 
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Figure 18: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Unspecific configuration is shown, pd data are from 
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Figure 19: nd break up cross section in [mb MeV"^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Unspecific configuration is shown, pd data are from 
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Figure 20: nd break up cross section in [mb MeV^^ sr^^] and nucleon analyzing 
power along the kinematical locus S (in MeV) at 65 MeV in comparison to predic- 
tions at NLO (light shaded band) and NNLO (dark shaded band) in chiral effective 
field theory. Unspecific configuration is shown, pd data are from 
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